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This paper extends the results of two previous papers on partitioned matrices. General 
reduction formulas are given for partitioned matrices A of order np satisfying A(XXI P ) 
= (XXI P )B, where B is a matrix of order rp and XXI V represents the direct product of an 
nXr matrix X of rank r with the identity matrix of order p. These formulas are related 
to the formulas given in the previous papers for partitioned matrices satisfying A(I P XX) 
= {I P XX)B. 



1. Introduction 

In a recent paper L. S. Goddard and H. Schneider 
[4] l discussed the relationship between matrices A 
and B, of orders n and /// respectively, which satisfy 



.LV XH 



(1) 



for some nXm matrix, X, of rank r>0. 

The results of Goddard and Schneider were gen- 
eralized in a paper by the author [5] for partitioned 
matrices A and B, in which corresponding sub- 
matrices (or blocks), A tj and B tji satisfy an analogous 
relationship to that given by (1); i.e., 



A ij X j =X i B ij , (i,j= 1 , 



,P) 



(2) 



where A fj is UiXn h X t is riiXnii of rank r u and B tj 
is rrhiXnij. 

If, in (2), each A tj is square, of order n, By is 
square of order r, and all X f are equal to an nXr 
matrix, X, of rank r, then A is of order pn, B of 
order pr, and (2) can be written 



A(I P XX) = (I P XX)B, 



(3) 



where (I P XX) is a pnXpr matrix which is the direct 
product of the identity matrix of order p with the 
matrix X; i.e., 



rx 



LXX-- 



X 



X 
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(The general definition for direct products is given in 
section 3.) 

It w r as shown by the author in [6] that, given a 
matrix A satisfying (3) for some matrices B and X, 
there exists a matrix /\ which can be constructed 
from X, such that 



'B 

P~ 1 AP= 



(4) 



so that pr of the roots of A are the roots of B. Thus, 
if r<Cn, or of B is reduced, (4) gives a reduction for- 
mula for A. Naturally a matrix satisfying (1) 
satisfies (3) with p=l. 

In theorem 2 of this paper a similar reduction 
formula to those for matrices satisfying (1) and (3) 
is given for partitioned matrices A of order pn 
satisfying 



A(XXI P )=(XXI P )B, 



(5) 



for some nXr matrix X of rank r, and the corre- 
sponding partitioned matrix B of order pr. The 
relationship between matrices satisfying (3) and 
(5) is discussed more fully in the next section. 

Other reduction formulas for partitioned matrices 
are contained in the results of B. Friedman [3], 
S. N. Afriat [1], and J. Williamson [10]. These 
formulas will be shown to be special cases of theorem 
2. 

The reduction formulas given in section 5 are 
applied in section 7 to certain special partitioned 
matrices. In particular, a 9X9 matrix, involving 
an arbitrary parameter jS, which arose in a problem 
on lattice points at the Bureau of Standards, is 
reduced by these formulas to a set of matrices of 
order 4, 2, 2, and 1 and from these reduced matrices 
all eigenvalues are found as functions of p. 



2. Inner- and Outer-Related Partitioned 
Matrices 

The formulas (1), (3), and (5) may be related by 
the matrix X. For instance, if the matrix A in (1) 
is a circulant matrix, 



a a\ 



&n-l a 



KJll 



&n-2 &n-l 



&n-3 a n-2 



tt n -l ^0 J 



(6) 



and X is the matrix of its characteristic vectors, 
then the related matrix satisfying (3) will have each 
of its submatrices a circulant of order n, and the 
related matrix satisfying (5) will be in circulant 
form; i.e., 



A= 



A.q A\ 



A. 



A n -\ A Q . . . A n 



^Ai A 2 . . . A 



(7) 



A matrix satisfying both (3) and (5) would be one 
which has the form (7), where each submatrix Ai 
has the form (6). D. E. Rutherford [9] has given 
reduction formulas for matrices of this type as well 
as for matrices related to tridiagonal matrices. 

We will call matrices satisfying (3) "inner-related" 
partitioned matrices, or i.r.p. matrices, those satis- 
fying (5) "outer-related," or o.r.p. matrices, and 
matrices satisfying both conditions "doubly-related," 
or d.r.p. matrices. 

We shall show that dual reduction formulas exist 
for o.r.p. and i.r.p. matrices and indeed that an o.r.p. 
matrix can be put into i.r.p. form by a simultaneous 
permutation of rows and columns. 

Since the formulas (3) and (5) are stated in terms 
of direct products, the general definition and three 
of the properties of direct products are given in 
section 3. 



3. Direct Products of Matrices 

W. E. Roth [8] derives some interesting results on 
direct product matrices and gives the general defini- 
tion for rectangular matrices, together with some of 
the properties. The three properties needed for our 
results are given after the definition below. Most of 
the properties of direct products can also be found 
in MacDuffee [7], although he deals only with square 
matrices. 
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The direct product of an mXn matrix, A, and a 
pXq matrix, B, in that order, is an mpXnq> matrix 
defined as follows: 



AXB-- 



Three properties of direct products: 

1. If the products AC and BD exist, 

(AXB)(CXD)=ACXBD. 

2. If A and B are nonsingular, 

(AXBy^iA-'XB- 1 ). 

3. If A is an mXn matrix and B is a pXq matrix, 

AXB=P(BXA)Q, where P and Q are obtained 
from the identity matrices I mp and I nq by 
simultaneous permutation of rows and columns. 



4. Similarity Between O.R.P. and I.R.P. 
Matrices 

Theorem 1: If A is an outer-related partitioned 
matrix of order pn, with square blocks of order p, 
satisfying (5), then A can be transformed by simul- 
taneous permutation of ?'ows and ^columns into an 
inner-related partitioned matrix, A, with blocks of 
order n, satisfying (3). 

Proof. Suppose A satisfies (5). Using property 
3 of direct products, given above, there exist non- 
singular permutation matrices, P and Q, such that 



XXI P =P(I P XX)Q. 
Then, from (8) and (5), 

A[P(I P XX)Q]=[P(I P XX)Q]B, 



(8) 



or 



(P-iAP)(I P XX) = (I p XX)(QBQ-*). 
So, if we let 



A=P~ l AP, 



B=QBQ~ 



the relationship (3) holds for A and B, and, since P 
is a permutation matrix, A is obtained from A by 
simultaneous permutation of rows and columns. 
From this theorem w^e can reduce any o.r.p. 
matrix, if we know the reduction formula for the 
i.r.p. matrix, by first permuting the rows and 
columns and then using the i.r.p. reduction formula. 
However, it is usually much simpler to use the 
direct reduction formula for o.r.p. matrices which is 
given in the next section. 



5. A General Reduction Formula 

Theorem 2: Suppose we have a matrix A of order 
n satisfying (7), or a matrix oj order pn satisfying 
(8) or (5) for some corresponding matrix B, where X 
is a given nXr matrix of rank r. Then in each case 
/lure exists a matrix P, depending upon X, such that 
(/+) holds, so that B contains r (or pr) of the roots of A. 

The only new result here is the formula for a 
inn (fix satisfying (5), as the formula for a matrix 
satisfying (1) was given by Goddard and Schneider 
[4], and the one for a matrix satisfying (3) was given 
by the author in [6]. The three transformation 
matrices are given below, however, in order to show 
the relations among them, but the proof is not given 
in detail as it appears in the other papers. 

Proof. Since X is an nXr matrix of rank r, 
we will assume 

x -Q m 

where Xi is a nonsingular (rXr) matrix, for if a 
permutation of the rows of X is necessary to achieve 
this, it can be matched by a simultaneous permuta- 
tion of the rows and columns of A] i.e., a permutation 
of the rows of X is equivalent to left multiplication 
by a nonsingular matrix V- but a matrix A which 
satisfies (1) also satisfies 



QAQ-\QX) = (QX)B. 



(10) 



In each case below we will use tin 4 matrix R, 
which we construct as follows: The first r columns 
of R are those of the nXr matrix X } which we assume 
to have the form (9). The remaining n—r columns 
consist of an rX {n—r) block of zeros and the identity 
matrix of order n — r, i.e. 



R- 



X\ 
Z 2 







h-r 



(11) 



Then 



R~ 










— X 2 X 1 1 I n _ r 



(12) 



Using R we can construct the transformation 
matrices for nonpartitioned matrices satisfying (1), 
or partitioned matrices satisfying (3) or (5) or both. 

Case 1 . Non-partitioned matrices 

Suppose A is an nXn matrix satisfying (1) where 
X is an nXr matrix of rank r and B is an rXr 
matrix. It was shown by the author in [6] that (4) 
holds for the matrix A if we let P=R. 

Case 2. I.R.P. matrices 

If A is a partitioned matrix of order np, with 
blocks of order n, satisfying (3), A also satisfies (1) 
where the matrix I v X X replaces the matrix X. In 



this case it was shown by the author in [6] that, if 
we let P=I P XR, then 

P~ 1 AP=A, 

where A can be put into the reduced form on the 
right side of (4) by a simultaneous permutation of 
rows and columns. In this case B is of order pr. 

Case 3. O.R.P. matrices 

If A is a partitioned matrix of order np, with 
blocks of order p, satisfying (5), then A also satisfies 
(1) with X replaced by 



(xxi P y- 



x x xi P 
x 2 xi p 



where, by property 2 of direct products, since X t is 
nonsingular, XiX/j, is nonsingular. Thus, if we let 



P=BXL 



XiXlj, 

\2Xlp l n - r XJ-i 



(13) 



(4) will hold in this case also, and B will again be 
of order pr. 

Case Jj. D.R.P. matrices 

If A is a matrix of order np, with blocks of order 
Uj and satisfies both (3) and (5), i.e. 



and 



A(XXln) = (XXl n )B, 

A(I P XY) = (I P XY)C, 



\\ here X is pXr and Fis 7iXs, then A can be reduced 
twice. For, by case 3, if we let P 1 = R i Xl n , where 
Ri is the pXp matrix given in case 3, we have 



P^AP,= 



II 



n 



where B is of order nr and D is of order n(p>—r). 

It is clear that premultiplication and/or post- 
multiplication of a partitioned matrix A , with blocks 
of order n, by a matrix whose blocks are scalar 
multiples of the identity matrix of order n produces 
a matrix whose blocks are linear combinations of 
the blocks of A. 

Also, if any set of matrices has the property that 



then 



or 



AtX^XBi, (1=1, . . . ,t) 

J2(c i A i X)=J2 ( Ci X B t ), 



(I] c € Ad X=X(J2c i B<), 



so that any linear combination of the set satisfies 
an equation of the form (1) with the same matrix X. 



57 ( j(M;0 — 61- 



Hence the matrices B and D satisfy 

B(I r X7) = (J r XY)E, 

D(I p _ r XY) = (I p _ r XY)G, 

where E is of order rs and G is of order s(p—r). 
Then, by case 2, there exist matrices P B and P D such 
that 

IB 

P^BP B =( 
\0 



PB l DP D - 



G 





2? 



Thus, the roots of .A are those of the matrices E, F, 
G, and #, of order rs, r(n—s), s(p—r), and {n— s) 
'(p—r), respectively. 

It would be easier in this case however, to multiply 
the two transformation matrices R } Xl n and I V XR 2} 
where R 2 is the nXn transformation matrix corre- 
sponding to the matrix in case 2, and obtain, by 
property 1 of direct products, the single transforma- 
tion matrix, P=R 1 XR 2) which could be used to put 
A into the reduced form, 



P~ 1 AP= 



rj^ # * * ^ 

F * * 

G * 

,0 5; 



For example suppose 



S= 



c 






where A, B y C, and D are 3X3 stochastic matrices, 
i.e. 

3 












(i=l,2,3) 



then S is an i.r.p. stochastic matrix as discussed in 
reference [6]. 
If we have also 

A+B=C+D, 

then aS is a d.r.p. stochastic matrix. 



Since a stochastic matrix with row-sum s has one 
root s with corresponding eigenvector (1, 1, . . ., 1) 
(cf. A. Brauer, [2]), the matrices R x and R 2 would be 



Rr- 



Then 



and 



P=R 1 XR 2 



p-^Rr'xR, 1 -- 



R 2 = 

R 2 
R 2 



r l 





°"1 


1 


i 





u 





u 





R 2 



R; 1 



So 



P-'SP^ 



R 2 1 (A+B)R 2 




where 



fa+b 




#22 #12 
#32 #12 



^2+^2 






R 2 l BR 2 
R 2 \D-B)R 2 

b 

5 2 

D 2 -B 2 J 



#23 #13 
#33 #13 



and B 2 and Z> 2 have similar forms. 

Another example using i.r.p. and o.r.p. transfor- 
mation matrices is given in section 7. 

In each of the above cases, if r=n, the matrix P is 
much simpler; i.e., for 1, P=X, for 2, P=I P XX, 
for 3, P=XXI P and for 4 if s=p, P=XXY. 



6. Theorems by Williamson, Afriat, and 
Friedman 

It was shown by the author in [6] that Williamson's 
theorem [10] on partitioned matrices follows from 
theorem 2, part 2, (theorem 1 in [6]), where r=n 
and the matrices B tj are triangular. 

S. N. Afriat [1] gives a general discussion with a 
number of theorems on the determinant and charac- 
teristic roots of partitioned matrices, A=(A ij ), in 
which the submatrices A tj all commute. One of his 
results is a generalization of Williamson's theorem 
and shows that if the roots of the commutative 
matrices A tj are \\f (k=l, . . .,n; ij=l, . . .,p) then 
the roots of the partitioned matrix A= (A i3 ) are 
roots of the matrices X A; =(X f ( - :) ). Since, as Afriat 
proves, commutative matrices can be simultaneously 
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reduced to triangular form by the same unitary 
transformation, it is clear that A=(A tj ) satisfies (3), 
where X is the unitary matrix in question, and the 
corresponding matrices, B ijf are triangular of order n. 
So Afriat's result on the characteristic roots of A also 
follows from theorem 2. 

Theorem 2 also provides a new proof for the 
following theorems by Friedman [3]. 

Theorem I. Let A u A 2 , . . ., A t be n-dimensional 
square matrices and B u B 2y . . ., B t be p -dimensional 
square matrices. Suppose that A u A 2 , . . ., A t have a 
common eigenvector, f , and the corresponding eigen- 
values are \ l9 A 2 > . . ., \ t respectively. Then p oj the 
eigenvalues oj 

C f =B 1 X^i+5 2 Xil 2 + . . . +B t XA t (14) 

will be eigenvalues oj the matrix 

D=\ l B l + \ 2 B 2 -\- . . . + \ t B t . 

As Friedman points out, this is a generalization of 
the well-known theorem that the eigenvalues of 
BXA are the products, X«/ij (i=l, • • *, n \ i =i ? 
. . .,p), where X* are the roots of A and ^ are the 
roots of B. 

Theorem II. Let C be the matrix considered in 
theorem 1 . Suppose the ring generated by the matrices 
A u . . .,A t has an r -dimensional representation (r<Cji) 
in which the matrix A k is represented by M k . Then pr 
oj the eigenvalues oj C are eigenvalues oj the matrix 

D=B 1 XM 1 +B 2 XM 2 + . . . +B t XM t . (15) 

Since theorem II contains theorem I as a special 
case with r=l, we will only prove theorem 1 1. 

Proof. Friedman shows in his proof that the 
statement, "A k is represented by M k " is equivalent 
to 



A k X=XM k 



(*=1, 



.,*) 



(16) 



where the columns of the matrix X are the basis 
vectors for the subspace. 

Let the elements of the matrix B k be b\f (i,j=l, 
. . .,p;fc=l, . . .,*). _ 

Then from (14), C is a partitioned matrix of order 
vn with blocks of order n, 



k = l 



(ij=l, . . -,p). 



(17) 



From (16) 
i±b<®(A k X)=^V${XM k ) 

k=l fc=l 



=X T.bfjM, (i,j=l,...,p). (IS) 

/c=l 



So if we let 



t-i 



(i,j=l ,p) 



(19) 



then D=(Dij) is the matrix defined by (15) and we 
have, from (17), (18), and (19) 



C i) X=XD il (i,j=l,..., P ), 



(20) 



so that the matrices C and D satisfy (3). Thus, by 
part 2 of theorem 2, there exists a matrix P such that 



-*m?;> 



(21) 



so that pr of the roots of C are roots of D. More- 
over, if the matrix X is known, D and E can be found 
explicitly by using (11) to construct P. 

7. An Application to a Lattice Point Problem 

The example below is a matrix which arose in a 
problem involving lattice points at the National 
Bureau of Standards: 



sf= 



4 -1 


-1 


£ 








-\ 


-1 4 -1 


~P 


-1 


P 








. -1 4 




~P 


— 1 


■ . 


. 
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-1 -0 


4 


-1 




-1 


P 


. 


(3 -1 -j8 


-1 


4 


-1 


~P 


-1 





. -1 




-1 


4 




~P 


— 1 




-1 


~P 




4 


-1 






£ 


-1 


~P 


-1 


4 


-1 






P 


-I 




-1 


4. 


we partition . y 


into 


3X3 


submatrices 


we see 



If 

(A B 0\ 
./=[B T A BY (22) 

\() B T A/ 

D. E. Rutherford [9] shows that a tridiagonal 
matrix, 

Ca b "N 

6 a b 

b a b 



T= 



b a b 
b a j 

has characteristic roots which may be written in the 
form \ k =a-{-2b cos 8 k where d k =kT/n J rl ) k=l, . . .'n. 

The characteristic vector corresponding to \ ; - is 

Zfc=(sin 6 k , sin 20 k) . . ., sin n6 k ). (23) 
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The transformation for an i.r.p. tridiagonal matrix 
was given by the author in [6] and Rutherford also 
gives all the roots in the d.r.p. case. 

The example (22) falls into none of the above 
categories, however, since B^B T ; but we show how 
the transformations discussed in section 4 can be 
used to simplify the matrix so that five of the roots 
may be found without having to solve an equation 
of degree higher than 2. I am grateful to D. E. 
Rutherford of the University of St. Andrews for 
giving me the other four roots. 

We first use the transformation matrix, Pi=(XX D 
where X is the matrix of vectors (23) for a 3X3 
tridiagonal matrix; i.e., 



X=- 



V2 

1 



V2 l" 

-V2 
■V2 1 



x- i =x, 



and / is the identity matrix of order 3 . Then we find 
r A-^2l C T 

^=Pf 1 j/Pi= C A C T 

^ D C A+J2IJ 



where 



C-- 



B-B T 



( p (T* 
-|8 p 
-0 OJ 



Since A and also the matrices A±-^2 I are all 
tridiagonal, we now use the transformation matrix 
P 2 = (IXX) and find 



<tf=P; l BP 2 = 



J9-V2 I C T 

C D C T 







O D+V2 IJ 



where 



^=diag (4-V2, 4, 4+V2). 



If we write the rows and columns of <& in the order, 
826435791, we obtain the matrix 



D= 



E 



F. 



with 



E= 



6 H 
HG 



F= 



K L 



where 



G- 



K= 



4+V2 p 
P 4-V2 

r 4 -0 . ■ 

-p 4 -p 

.-pi. 



H= 



r 



(T| 
d, 



and M=diag(4+2-v/2, 4 — 2 A /2). Thus £ is an o.r.p. 
circulant, so the roots of E are those of G-\-H and 
^— fl", which are 4±/3±ViS 2 + 2. 
Since i£is tridiagonal, let 



/X 



ft= 



Then 



P 3 -^3 = 



■4-V2/? 



/3/V2 



^/V2 



/3/V2' 



4+V2/3 -0/V2 -/5/V2 
-18/V2 4+2^2 . 
-18/V2 . 4-2^/2^ 



Thus 4 is a root of this matrix, and when we 
delete the second row and column, we are left with 
a 4X4 matrix. I was unable to reduce this further 
by any relatively simple similarity transformation. 
However, D. E. Rutherford has found the other 
four roots of F: 



4±V/3 2 +2/H-2±V/ 52 -2/3+2. 
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